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Abstract 

This paper considers some physically interesting cosmological dynamical 

systems in the FRW-scalarfield category which are examined for integra- 

^H, bility according to the criterion of Painleve. In the literature these systems 

have been examined from the point of view of dynamical systems and the 
results from the two disparate methods of analysis are compared. This 
allows some more general comments to be made on the use of the Painleve 
method in covariant systems. 
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c3 ■ 1 Introduction 



This is the seventh in a series of papers ||i|, ||, y, H B bI devoted to various 
aspects of the three mathematically disparate topics of symmetry, singularities 
and integrability of dynamical systems. It is the purpose of the series to cement 
the concepts into a greater obvious unity whereby they are recognised as dif- 
ferent aspects of the central properties of dynamical systems vifhich distinguish 
integrable systems from nonintegrable systems. 

The plan of the present paper is as follows. In §2 we present a brief resume 
of the analysis of Painleve with particular reference to some of the subtleties 
required to deal with systems which do not fit into the neat scheme of the 
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standard ARS algorithm |Q ^, ||, ^. In §§3-6 we present the analysis of a 
number of models chosen from some interesting areas of Cosmology which have 
been analysed from the viewpoint of Dynamical Systems. In the approach used 
in Dynamical Systems the original system of equations is replaced by an n- 
dimensional system of first order equations, maybe subject to a constraint to 
preserve the actual number of degrees of freedom. This is not always a suitable 
format for the application of the Painleve analysis and we find it necessary to re- 
arrange the systems into a more suitable format. The possession of the Painleve 
Property, i.e., the existence of a solution analytic except at moveable polelike 
singularities (branch point singularities in the case of the 'weak' property), is 
representation dependent, being preserved under a homeographic transforma- 
tion only, and consequently we must be careful in comparing the results obtained 
using the Painleve approach with those obtained using the methods of Dynam- 
ical Systems. Finally, in §0, we present some general observations of the results 
of our analysis paying particular attention to the uses of Painleve methods and 
analyses when dealing with general relativistic systems. 

2 Methodology of singularity analysis 

For a detailed account of our approach to the Painleve analysis the reader is 
referred to our earlier paper B . Here we present a summary of that account so 
that the paper is self-contained for the reader not concerned with the fine detail 
of the ideas behind the singularity analysis. 

The basis of the singularity analysis of a system of differential equations 
is found in the dominant behaviour of the system at a singularity and also 
the next to leading order behaviour in the neighbourhood of the singularity 
[pi O. When the analysis is concerned with the integrability of a system, 
the conditions on the behaviour must be tightened to analytic behaviour in the 
complex plane or at least substantial portions of it. (One can treat with branch 
point singularities provided the branch cuts do not become too concentrated in 
the region of the complex plane relevant to the solution of the problem at hand.) 
This was the original approach by the originators of this field - Kowalevskaya 
[ [l3| and Painleve ||lj, |l5[ |6|. The idea behind the analysis is that a system 
of differential equations is integrable if there exists a set of analytic functions 
which satisfy the system of differential equations. Since an analytic function 
can be expressed in terms of Laurent series, the demonstration that there exist 
Laurent expansions about movable singularities containing the requisite number 
of arbitrary constants for the dependent variables in the system of differential 
equations is taken to indicate that the system is integrable. This is the basis 
of the ARS algorithm developed by Ablowitz, Ramani and Segur [0, |[ ^ and 
presented, for example, by Ramani et al [jl^l and Tabor [ p7[ . 

The essence of the ARS algorithm is that the solution of an n-dimensional 
system of equations 

x,^f,{t,x), (2.1) 



where the functions fi are rational in the dependent variables and algebraic in 
the independent variable, can be written as either 

x,(i)=^a,r-P*+^"' (2.2) 

j=o 

in the case of a Right Painleve Series or 

oo 

x,(i) = ^a,r-^-^-« (2.3) 



in the case of a Left Painleve Series [^ 11 , where t ~ t — to and to is the 
location of the movable singularity, the exponents pi are strictly positive integers 
(rational numbers in the case of the so-called weak Painleve test), the paramater 
g is a positive integer (respectively rational number) and in the coefficients there 
are n— 1 arbitrary constants which, together with to, give the required number 
of arbitrary constants for the general solution of the system (^J). In the case 
of the Left Painleve Series, which has the nature of an asymptotic expansion in 



the region of t infinite (p.3[) is usually written as 



xdt)^ya,t-P'-^^ (2.4) 






and the coefficients now contain n arbitrary constants. 

The ARS algorithm provides a mechanistic approach to the establishment 



of the existence of (2^) and/or (p.4|), which is not always satisfactory ||19|, but, 
when handled with sensitivity rather than automatically, provides a useful tool 
in the analysis of differential equations. The first step is to determine the leading 
order behaviour by means of the Ansatz 

x^^a^T-P% (2.5) 

which is substituted into the system (^) , and to assemble all possible patterns 
of singular behaviour. For each of these possible patterns the coefficients at are 
calculated. For each pattern the substitution 

X, = a,T-P^ + /i,r'-P- (2.6) 

is made to determine the 'resonances' r at which the arbitrary constants, fii, 
are introduced. These are determined from the leading order terms of (p3|) by 
a linearisation process. 

For each pattern of singularity behaviour and its concomitant resonances 
and arbitrary constants the substitution 






a.r^P'+^«, (2.7) 

3=0 



where jmaxl = fmax and Vmax is the largest resonance, is made into the full 



system (2.1) to ensure that there is consistency among the resonances and with 
the nondominant terms of the original system. If this be the case, that pattern 
of singular behaviour passes the Painleve test and, if this be the case for all pos- 
sible patterns of singularity behaviour, the system (p^) possesses the Painleve 
Property which guarantees integrability in most cases. There are exceptions 
to this general statement and the interested reader is referred to Miritzis et al 
[H and the references cited therein for a more detailed discussion of the finer 
points of the application of the ARS algorithm to the techniques of singularity 
analysis. 

There are systems which cannot possess the Painleve Property since they 
do not have movable singularities in every variable. We speak not of those for 
which one or more of the variables has a leading order behaviour in terms of a 
positive rational exponent. This becomes singular after a sufficient number of 
differentiations. Rather we refer to those variables which are constant at the 
movable singularity, toj of other variables. (Positive integral behaviour can be 
rectified by the simple, yet homeographic, means of taking the inverse as a new 
variable.) Usually an analysis of these cases by means of the Ansatz of a series 
expansion - the ARS algorithm in not applicable in such a case and, indeed, is 
downright misleading ~ produces series with fewer that the required number of 
arbitrary constants. Naturally such series cannot represent the general solution 
and give no indication of the integrability of the system. However, subject to 
the convergence of the series, they can indicate some sort of partial solution. For 
reasons explained in Miritzis et al M we refer to these solutions as 'peculiar' 
solutions - a literal translation from the Greek usage - to avoid some of the 
associations with other terms commonly used such as 'singular' and 'particular'. 



A peculiar solution is one consisting of series which satisfy the system (2.1) 
containing fewer than the required number of arbitrary constants. Finally we 
note that t he co nsiderations presented here in terms of the systems of first order 



equations, (2.1), apply mutatis mutandis to all systems of ordinary differential 
equations. 

3 Massive scalar field 

Consider the Friedmann equations for the evolution of the scale factor a{t) of 
an open, flat or closed {k = —1,0,-1-1 respectively) universe in the case where 
matter is represented by a self- interacting scalar field with potential T^ (</>). They 
are 



(3.1) 



\aj a^ 


3 V2 J 


2^ = 
a 


= -1(20^-^(0) 


3-0 + r(0) = 

a 


= 



In this Section, we focus in the particular case that of the scalar potential, V {(/)), 
is quadratic of the form ^nv^cjp so that the above equations become 



1 /l :^ 



a 1 



—m 
2 



6 



202 



2 J.2 



-m (f> 



(3.2) 



+ 3-(/> + 771^0=0. 

a 



These three equations are not independent since, for example, the first and 
third of (3.2) lead to the second. (This is also the case in (3.1). The result is 
independent of the precise form of the scalar field.) We make the choice (3.2c) 
and so treat the system 



k 



m^cj)^ 



3-0 + m^0 = 0. 
a 



(3.3) 



We can make some 'cosmetic' improvements to (3.3) by rescaling time, t, scalar 
field, (j) and scale factor (in the case k ^ 0), a, according to 



t = mt. 



to obtain the system 



iiiv 



¥' 



u 
uv 



and u — 



ma 



u" {v^ + v^) 



in the case k = Q and the system 



Au = Au^ (v^ 



2uv + Siiv + 2uv = 



(3.4) 



(3.5) 



(3.6) 



in the case that k ^ 0. In both instances the overdot now represents differenti- 
ation with respect to the new time variable, i. We observe that the sign of k is 
immaterial in the transformatio n (t3.4| )- 

We commence the testing of ( 3.5) for the possession of the Painleve Property 
with the standard substitution 



ar 



Pt\ 



(3.7) 



where t — t — to and to represents the location of the moveable singularity if p 
is found to be negative, to determine the leadin g or der behaviour. We obtain 
the following exponents for the several terms in ( |3.5| a) and (3^b) respectively 



2p-2 2p + 2q-2 

p + q — 2 p + q-2 P + Q- 



2p + 2q 



(3.8) 



From (3.8a) we see that p is arbitrary and that balancing between the first 
and third terms is possible for q = —1. However, then the second term has 
the lower exponent 2p — A and cannot balance if one is thinking in terms of a 
Right Painleve Series. However, there is no problem if one thinks in terms of 
the existence of a Left Painleve Series since now the second term vanishes more 
rapidly than the two dominant terms at infinity as required. Unfortunately the 



third term of (3.8b) is incompatible with a Left Painleve Series and one must 
conclude that q = which thereby ends the validity of the next to leading order 
analysis. The standard procedure outlined in the ARS algorithm ceases to be 
applicable. 

For { \i.(\ ) we apply the same methodology as for (3^) and obtain the same 



dismal results. There is no possibility of consistent singular behaviour. 

In the cases of both ( ^.5[ ) and (3^) one can look for the existence of a 
standard series solution in the absence of the Laurent series implicit in the 
Painleve analysis. Indeed we find 



ao + ait + 



3albl 



2ao 



2a? ^ 



bo 



^/aj-apl - 



ao 



bo 



3 — Oi 



P 



(3.9) 



in the case of (3.5) and 

u = ao + ait + I 3aofoo 

A /Ann - 

V = bo 



Aa, 



P + 



V^ 



2l2 



4ao^o 



2ao 



t- 



](2bo + 3—h]P + 
A \ ao 



(3.10) 



in the case of (3^). The coefficients ao, bo and ai are arbitrary in both cases. 
Coefficients of higher powers of t become somewhat more complicated expres- 
sions. 

We conclude that t he model of the massive scalar field with a quadratic 
potential described by ( |3.3|) is not integrabl e in the sense of Painleve no matter 
the curvature (the change of variables in (3^) can at most have the effect of 
a change from possession of the Painleve Property to possession of the 'weak' 
Painleve Property) . This is of course not surprising given that very little is now 
about the evolution of this system in general. 

However, making the usual slow-roll approximation, which basically amounts 
to neglecting the (p and kinetic energy terms from the Friedman equations, that 
is 






(3.11) 



implies that in the quadratic potential case we are considering, the exponents 
may balance with p arbitrary and q = —1. This again docs not obviously 
produce some great difference in the possession of the Painleve property but 



it does show that at least the ARS algorithm can be applied. So we see that 
precisely this condition, which is conducive to inflationary solutions, makes some 
of the difficulties disappear. 



4 Models with an exponential potential 

4.1 The scalar- vacuum case 

The equations for the simplest case, that of an FRW model with a single scalar 
field with an exponential potential V = e"'^'* and no other matter source, are 
|0|: 

a = -d^ - 2(p^ + -Xa(f) + 1 



= iA(/.2-3d0 + iA 



subject to the constraint 



l-ke 



-2a+\cf> 



(4.1) 



(4.2) 



where, as usual, /c = 0, ±1 depending upon whether the scalar curvat ure is zero , 
posit ive or negative. As we noted in §|| in the instance of equations (3J), (4.1) 
and (O) are not independent sin ce the differentiation of(4.2) and one of (4.1) 
leads to the other. However, (4J) and one of (11) and (12) are not equivalent 
systems since one is of the fourth order and the other of the third order. 



4.1.1 The scalar- vacuum case as a third order system 



For our third order system we take (lib) and (12), i.e., we consider the system 



1 



J,2 



1 



X6^ ~3a(j)+-X 

2 2 

6^ + 1- fce^2a+A0^ 



(1.3) 



We transform the system (4.3) into one suitable for application of the singularity 
analysis by means of the nonhomcographic transformation 



-2a+A0 



oA0 



to obtain 



2uv — 3uv — X uv = 



■ 2 2 
U V 



2uuvv 



1 - — I uH^ = 4m^w^ 



Aku'^v 



3„,2 



The dominant expressions in the leading order analysis are 
2a[3q{q - l)rP+«~2 _ 2,aPpqTP+''-'^ 



(4.4) 

(4.5) 
(4.6) 

(4.7) 



and 

^2^2p2^2p+2,-2 



2a''l3^pqT 



2p+2g-2 



4 
A2 



^2^2^2^2p+2,-2 ^ 4fc^3^2p2^3p+2, 

(4.8) 



for (4.5a) and (4.5b) respectively after the usual substitutions are made. From 
(4.8) we obtain p =^ —2,5 arbitrary and, using this in (4.7), we find that q = 
-2,0. 



Firstly we consider the pattern of singular behaviour (—2, —2). From (4 
we obtain a = 4/fcA^ and /3 remains arbitrary. When we make the substitution 

-2 I 1 — 2 a —2 I r — 2 

the condition that the required constants are arbitrary is that 



(4.9) 



6r 
16 

"A2 



2r(r - 2) 
16r 

^2" 







(4.10) 



which has the solutions r = —1,0. The first is generic, the second reflects 
the arbitrariness of /3 and there is no third so that the system cannot pass the 
Painleve test with this pattern of singularity behaviour. The solution is peculiar 
(c/ ID) with the first few terms being given by 



(4.11) 



which is suggestive of the solution of (4.5) when v — constw and the pole is fixed 
at t = 0. 

The singularity pattern (—2, 0) is without the ambit of the Painleve test. If 
we make the Ansatz 




i=0 



aiT 



4=0 



(4.12) 



the requirement that the coefficients of the first two terms balance gives 



ao&i = 



2 1,2 
an On 



kal% = 0. 



(4.13) 



The choice oq = 0, 6o 7^ leads to nonzero even coefficients for the bi and for all 
Oi = 0. The reverse choice gives the only nonzero coefficient to be oq. Clearly 
both are just nonsense solutions. 

From the point of view of a third order system the scalar- vacuum case does 
not pass the Painleve test, does not possess the Painleve Property and so is 
not integrable in the sense of Painleve. However, we do note that there does 
exist a peculiar solution of the type mentioned by Ince pl^, p 355] in the context 
of an equation discussed by Chazy in 1912. Naturally this solution would be 
acceptable were the system of the second order. This leads to the approach 
discussed in the following subsection. 



4.1.2 The scalar-vacuum case as a second order system plus a first 
order system 



An alternate approach to the treatment of systems (4.1) and (|4.2|) is to regard 



(4.1) as the two dimensional system of first order equations 



-x' + -Xxy - 2?;2 + 1 

-3xy+iAy2 + iA, (4.14) 



where x = a and y = <f). This system can be treated independently of (4.2) 
since the latter contains the two constants of integration following from the 
integration of x and y. (We note that this approach was used by Halliwell p3 



to treat (4.14) as a two-dimensional dynamical system.) 



The usual substitution for the leading order behaviour of (4.14) gives the 
following patterns of exponents 

p-l 2p p + q 2q 
q-l p + q 2q. 

The values p = q = —1 are consistent for both sets. This is the only possible 
pattern since, although p = —1, q > — 1 is fine for the first set, it fails for the 



second since the left side of ( 4.15 1) becomes zero and so no balancing can occur 
and the reverse happens with p > —l,q = — 1. The coefficients of the leading 
order terms satisfy 

-a = -a^ + -XaP - 2/3^ 

-f3^~3a(3+^\/3^ (4.16) 

from the second of which it follows that 

« = i + iA/3 (4.17) 

and from the first of which with ( [4.17| ) 

(A^ - 36) /32 + 4A/3 + 4 = (4.18) 

We recall that A > 0. Consequently the coefficients of the leading order terms 
are given by 

a=i/3 = -i if A = 6 (4.19) 

6 6 

" = TA /3=T^ if 0<A^6. (4.20) 

A± 6 

This means that for general A there are two possible solutions to be considered 
for consistency with the requirements of the Painleve test. 



The resonances are the eigenvalues of the equation 



1 



-2a 
3/3 



P\ -iAa + 4/3 
r-l + 3a~ \p 



which foUows from the substitution 






= 



(4.21) 



fir 



and y — (3t ^ + , 



(4.22) 



into the dominant terms of the system ( 4.14 ). We obt ain th e generic ri = — 1 
no matter the value of A and r2 = | in the case that ( 4.19 ) applies and r2 — 
±8/(A±6) = t4/3 in the case that ( 4.20 ) applies. Since the nondominant terms 



in (4.14) are constants, the expansion must be in a power of t such that it be 
commensurate with 1 and the second resonance, r2. This involves a constraint 
on the permissible values which A may take. It is impossible for both values of 
r2 to belong to the interval (0, 1) and so in this range the system (4.14) cannot 
pass the Painleve test. However, the requirements of commensurability and two 
positive values of r2 can be satisfied for some A G (0, 6) and the potential of both 
branches to pass the test is present. The possession of the Painleve Property 
by (4.14) is not possible in the case r2 = 1, i.e., A = 2 since at that value the 
arbitrary constant that is introduced has zero coefficient and cannot be used to 



compensate for the presence of the nondominant constant term in (4.14). The 
constant terms deny the possibility of a Left Painleve series and are responsible 
for the restriction on the interval of A for which interval one can hope to obtain 
Right Painleve series for both solutions. 



Note that possession of the Painleve Property by (4.14) means that the 
original second order pair ( [4.l| ) does not possess it since both a and cj) will 
contain a logarithmic singularity due to the quadrature of the leading order 
terms. Nevertheless a solution of lessened an alyt ic quality will exist and one 
must now consider the e ffect of the constraint ( [4.2[ ) . After a little manipulation 
the combination a:( 4.14 a) — y(4T4b) leads to 



xi-yy=- {x^ - y^ - l) {-2x + Xy) 



which is immediately integrable to 



y 



1 — est exp(— 2a + X(f>). 



(4.23) 



(4.24) 



Consequently there is a further constraint on the coefficients of th e exp ansions 
for x{t) and y(r) determined by the constant of integration in ( 4.24 ). This 
constant is the value of fc, = 0, ±1, and this will fix the single arbitrary constant 
found in the series (apart from the location of the singularity) . 

For example, in the case that A = 10 so that r2 = 1/2, the expansions are 



x{t) 



ri=0 



n/2-1 



yW 



n=0 



(4.25) 



10 



The result is surprisingly nice as the first few coefRcients are 



ao = 1/8, 


bo = -1/8 




ai = 36i, 


bi^bi 




02 = 26?, 


62 = -26? 




03 = -246?, 


bs = -86? 




04 = i (11 - 1926?) , 


64 = i (13 +1926?) 


(4.26) 


and we obtain 






x(t) = i(t-to)-V36i(t-io)"'/' 


+ 26?-246?(i-io)'^' + 0(t- 


-to) 


yit)^—^{t-tor^ + b,{t-toy^^'' 


'-26?-86?(i-io)'/' + 0(t- 


-to). 
(4.27) 



Substitution in the last constraint equation seems to be a formidable task. The 
left hand side is 

LHS = 61 (t - to)'^^^ + 86? (t - to)"^ + 86? {t - to)~^^^ - 1 - 1286? 

- 1286? {t - tof'^ + 5126? (i -to) + 0{t~ tof'^ (4.28) 

and the right hand side, up to the constant of integration, 

RHS ^{t- to)'^^^ + 861 {t - to)'^ + 86? (i - to)^^^^ - 1286? 

- 4806? {t - io)'/' + ^bl (t -to) + 0{t- io)'/' . (4.29) 

Consequently and therefore we can only conclude that, although there is 
agreement in the two series for the initial terms, this agreement fails for higher 
terms in the Laurent expansions. In the neighbourhood of the singularity the 
two would appear to be almost equal and in a numerical integration one could be 
tempted to ascribe the discrepancies to numerical error and to input integrability 
in the sense of Painleve to a systems which is not integrable in that sense. (One 
recalls a similar suggestion of integrability in the surface of section plots for 
the Henon-Heiles problem at low energies, the equivalent to being close to the 
singularity in the present instance.) 

4.2 The scalar-fluid case 

This is the case where a general (flat or curved) FRW model has a scalar field 
(with an exponential potential) and a separately conserved perfect fluid. The 
equations for the flat case are |g^ [Eqns (6)- (7)] 

x' = -^(2 - 7)x + X^^y^ - 3(2 - j)x^ - 3jxy^ (4.30) 

y' = Ay - ^Jl^y - 3(2 - i)^''y - ^ly' (4.3i) 



11 



which are derived from the system 



subject to 



P1 = 


-SH{p^+Py) 


4> = 


-3iJ(/)-y'((/)) 


Pj^ 


(7 - 1)P7 


V ^ 


Voe"^'"^ 


H^ = 


1 2^ 1; 



V 



(4.32) 

(4.33) 
(4.34) 

(4.35) 



which represent the assumed equation of state, the assumed form of the po- 
tential and the constraint following from the Einstein Field Equations, by the 
definitions 



Ven' 



y 






In a 



(4.36) 



and the prime denotes differentiation with respect to new time, r. In terms of 
the solutions of the system (4.30) and ( 4.31| ) the constraint equation ( 4.35| ) now 
becomes the definition of the variable nJ^j^fiH and so the treatment of the 
system ( 4.30D and (4-31) can be undertaken without the necessity to consider 
any constraint. 
We substitute 



to obtain the set of exponents 



2/-/3r« 



p~ 1 p 2q 3p P + 2q 

g~l q p + q 2p + q 3q 



(4.37) 



(4.38) 



from which it is evident that there are two possible patterns of leading order 



behavior, viz p : 
equations and p 



0, q 



involving the left side and final two terms of both 
— -i involving the second and fourth terms of the 



right side of ( f.3C ) and the left side and final term of the right side of ( 4.31 ) . 
The second pattern does not fit into the scheme of the Painleve test and so the 
possession of the Painleve Property depends only upon the first pattern. We 
note that the third order terms in both equations can also balance since they 
share the symmetry xdx + ydy |T^ , but this leads to some constant solutions 
which are not comparable with the remaining terms in the two equations. 

Case p = q = —^^■ 

From the dominant terms of both equations the coefficients are required to 
satisfy the single equation 

7/32 = (2 - a)2 - 1, (4.39) 



12 



which is very encouraging for the passing of the Painleve test. To determine the 
resonances we substitute 



ar 



^IT 



_i _i 

y ^ (3t 2 + vt"" 2 



into the dominant terms to obtain the two dimensional system 

= 



r + l-37/32, 
-3(2-7), 



3a/37 
r + 37/32 



V 



(4.40) 



(4.41) 



in which (4.39) has been used to make some simplification. 

The system ( [4.41 ) has a nontrivial solution if r = —1,0 which is consistent 
with (4.39) and means that the system ( 4.30| ) and (4.31) possesses the (weak) 
Painleve Property since the nondominant terms do not enter until higher powers. 

For the other pattern of behaviour we substitute 



r2' 



and 



°° 1 1 



(4.42) 



z=0 



We find that the odd coefficients are all identically zero and that the only arbi- 
trary constant is the location of the movable singularity (in y), t^. Consequently 
we have a peculiar solution. 

In terms of the definitions of the new variables in terms of the original 
variables the constraint equation, ( 4.35| ), gives kJ^J^ZH as a function with 
a square root singularity at ig analytic, away from the single branch cut, in 
the variable r. Unfortunately the definition of new time makes all functions 
nonanalytic in terms of the original time variable. However, there is a sense 
of integrability provided one keeps to the new time. It is not possible to make 
even this remark about the original dependent variables since it is the three 
combinations of them which have been shown to have the moderately good 
behaviour. 

The equations for the case where a general curved FRW model has a scalar 
field with an exponential potential and a separately conserved perfect fiuid are 
given by van den Hoogen et al p3 [Eqns (1.10-1.12)]. From the fourth order 
system 

1 



H = --[jp^ + <!>'] -K 

P-y = —3jHpj 

4>= -iH4) + KV 



subject to the constraint 



H^ 



P-t + T;r + V]+K 



(4.43) 



(4.44) 
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we obtain the third order system 

x' = -2a; + Ay^ + Bx?L + 2x^ - xy^ 
y' = y - Axy + ByQ + 2x^y - y^ 
n' = -2Bn + 2BD? + Ax'^Q. - 2y^n, 



(4.45) 
(4.46) 
(4.47) 



in which we have written A — ka/ | and _B = ^7 — Ito make the expressions 
look sHghtly simpler. 

Under the usual substitution for the leading order behaviour we find the 
following pattern of exponents 



p ~ 1 p 2q p + r 3p p + 2q 
q - I q p + q q + r 2p + q 3q 
r — \ r 2r 2p + r 2q + r. 



(4.48) 



It is evident that the general singularity pattern is p = q = —^,f = ^1- There 
are also other patterns of a more specialised nature. We list them in Tabic I. 



p 


q 


r 


Observation 


i 

9 


>-k 


> -1 


not P 


>-^ 


1 
2 


> -1 


not P 


>-i 


>-^ 


-1 


not P 


1 

! 
2 


1 
2 


> -1 


Pifr = -i 


>-^ 


-1 


not P 


>-^ 


1 

5! 


-1 


not P. 



Table I: Possible sundominant patterns of semisingular behaviour 
for the system 4.45 - 4.47. 



Only the fourth case is a candidate for the Painleve test as the other pos- 
sibilities do not conform to the singularity requirement of the function and/or 
its derivatives. 

In the general case the coefficients of the leading order terms satisfy the 
single constraint 



Aa^ - 2/3^ + 2B7 +1 = 
so that one would expect to find the resonances r = 



(4.49) 
-1,0(2). With the use of 



(4.49) the characteristic equation for the resonances is 

r - Aa^ 2a/3 -Ba 

-AaP r + 2/32 -B/3 
-8a7 4/^7 r - 257 



<^ r" 



r^ = 



(4.50) 
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(4.47) 



so that r = —1,0(2) as anticipated from (4.49) so that there is consistency. 
There are three arbitrary constants and the nondominant terms can cause no 
theoretical difficulties. For this pattern of singular behaviour (4.45) 
passes the Painleve test. 

The other candidate for the Painleve test, case four of Table I, with p 
r = — i produces the system 

1 ., , 

a = 2a'^ - af3 

2 

-i/3 = 2a2/3-/?3 
--7 = 40^7 -2/3^7 



(4.51) 



for the coefficients of the leading order terms. There is a contradiction between 



the first two equations of (4.51) on the one hand and the third equation on 
the other. Consequently this is not an admissible case for the system (4.45) - 
(4.47). The one admissible singularity pattern has passed the Painleve test and 
so we conclude that the system ( 4.45|) - (4.47) possesses the (weak) Painleve 
Prope rty. The constraint equation, ( [4.44 ), for the original fourth order sy stem . 
( 4.43 ), defines the fourth variable. The same comments as for the system ( 4.32 ) 
apply to this case also. 



5 The general flat case 

We use the notation of Foster |24] . He presents the dynamical system 



dt 

dK 3-2 

^ = "2^ 



(5.1) 



subject to the constraint 



K^^W{cj,) + ^4>^. 



(5.2) 



Foster introduces the new variables 



K' 



x—^T^ V ~ \ 7,~f7 ^^^d T = vogv (t) + tq 



2K 



(5.3) 



15 



where K is the trace of the extrinsic curvature and v{t) represents volume, 
V ^ a^ . Then the system becomes 

dx 



dr " 


9 


dy 
d^ " 


= -y 



x'V {(b) + y' 



dr 



and the constraint is now 



y^ + 3x^V{<j>) = 1. 



(5.4) 



(5.5) 



In order to make some reasonable progress with the Painleve analysis of 
this system we require some explicit functional form for the potential, V{(j))- 
(We note that in Foster's Dynamical Systems approach |24| such a constraint 
was not necessary. All that was required were some general properties of the 
potential. The approaches via Dynamical Systems and Singularity Analysis are 
complementary.) We make the standard Ansatz for a confining potential (c/ 
the transition from (3.1) to (3.2)), viz 



so that the system (p^) is now 



and the constraint is now 





F(0) = 


'- im202 




3 now 








dx 
d^^ 


2 

y X 






dy 
dr ~ 


-'^-\ 


j-x^m^^ + 


y 


d(t) 
d7^ 


\lly 






y' 


'+r" 


m^(j)^ = 1. 





(5.6) 



(5.7) 



(5.8) 



We make the usual substitution for the leading order behaviour to obtain 
aprP-^ = a(3^TP+^'i (5.9) 

l2 



^2^2^2^2v+r 



r-1 

ryrr = 



/3t« 



(5.10) 
(5.11) 



From (5.9) we have: q — —^ a is arbitrary, p arbitrary and 0^ — p. With 
this ( |5.11| ) gives r = i and 7 = 2^/2/3/3. In ( |5.10| ) the possibihty that p < -1 
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cannot occur as there would only be one term and so there can be no balance. 
If p > —1, the two dominant terms require that /?^ = — ^ and from above this 
means p = — i but a remains arbitrary, li p = —1, a^ = —l/Am?. 

The system ( |5.7D permits the two patterns of dominant behaviour, viz 
(— ^, — ^, ^) and (— 1, — ^, ^). For the first pattern of dominant behaviour the 
resonances are found to be r = —1,-^,0. The first is the generic value. The 
third reflects the arbitrary value of a. The second demonstrates that this pat- 
tern cannot pass the Painleve test since it requires a Left Painleve Series and 



this is inconsistent with the nondominant terms in the original system (5.7) 
which require a Right Painleve Series. The series expansion for this pattern of 
dominant behaviour contains only two arbitrary constants and consists of even 
terms only. There are two expansions depending upon whether a positive or 
negative root is taken. The first few terms of the series based on the positive 
root are given by 



2„2^^-2 



_1 1 . „ 2 

X = aoT 2 H — flo (1 + 2m a, 

i _i iV2 , o ON _i 

y^^T 2-^(1 + 2m2a2) r 2 

(/.= -^T2-^(l + 2m2a2)T-2 (5.12) 

For the second pattern of dominant behaviour we find that r = —1(3) which 
means that there is only one arbitrary constant, the location of the moveable 
singularity. The series expansion consists of even terms only and for the positive 
root the first few terms are 

i 1 Zi 155-243i 

2to 16r7?, 9216to 

_i 3i 155 1 

^ 16 4608 

, „ /2" 1 i [2 1 31 /2 5 

6 = 2\ -it2 a/-t2 H '\/-t2 .... (5.13) 

^ V3 16 V 3 3204 V 3 ^ ' 

Neither pattern of dominant behaviour passes the Painleve test and so the 
system (pj) does not posses the Painleve Property. 

6 Scaling solutions 

Billyard et al J25J discuss the stability of cosmological scaling solutions within 
the class of spatially homogeneous cosmological models with a perfect fluid 
subject to a linear equation of state and a scalar field with an exponential 
potential. This type of study was extended by van den Hoogen et al |E3| 
to Robertson- Walker spacetimes with nonzero curvature. Miritzis et al [^| 
have already shown that the corresponding problems in flat space possess the 
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Painleve Property. For the former problem Billyard et al reduced the governing 
equations to the three-dimensional system 



y|,y2 ^ 3^. [(^ _ 2) ^ _ 2^2] 



x = M^^y 



-jV [(7 ■ 



y = 3y-d -Kxy + -y [(7 - 2) 17 - 2y 

h = -3(7 - 2)n + sn [(7 - 2)n - 2y^] . 



(6.1) 



System (6.1) has a superficial resemblance to a cubic system, but this re- 
moved by the nonhomeographic transformation y^ = z which reduces (6.1) to 
the quadratic system 

^i^z + ^x[{^-2)n-2z] 
i = 6z — VGkxz + 3z [(7 — 2)il — 2z] 

n = -s{j-2)n + 3n[{j-2)n-2z]. (6.2) 

Since the transformation between the two systems is not homeographic, the 
possession of the Painleve Property by the one does not guarantee its possession 
by the other. However, possession of it by ( |6.2| ) will mean that {6A) will possess 
the so-called weak Painleve Property since the transformation replaces a polelike 
singularity in z with a branch point singularity in y. 
The standard substitution 



X = arP, z = (iT'i, Vt = Ft'' 



(6.3) 



into the derivatives and quadratic terms of (3.2) gives the set of exponents 



p— 1 q p + r p + q 
q — I P + q q + r 2q 
r — 1 r 2r q + r. 



(6.4) 



The third of ( |6.4| ) gives q — r = —1 and the second p = —I. The coefficients of 
the leading order terms are found from the solution of the linear system 



-(7- 
2^' 



2)aF - 3a/3 

-P = y/QnaP + 3(7 - 2)/3F 
-F = 3(7-2)F2-6/3F. 



6F^ 



(6.5) 



For nontrivial coefficients a, j3 and F we see immediately that there is an in- 
consistency between the first and third of ( |6.5| ). This can be resolved by setting 
either p = or r = 0. However, the latter choice leads to further inconsistencies 



and we are forced to the set p = Q, q = r = —1 and so the system (6.1) cannot 
possess the Painleve Property. 
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One could substitute the partially singular expansion 

x^J2''^^'' ^/ = 5]&.r(*-l)/^ r! = ^c,T'~i (6.6) 

i=0 i=0 1=0 

into ( |6.lD . The results are somewhat disappointing. The coefficients in the 
expansions can be expressed in terms of the one arbitrary coefficient, 69 • In 
particular 

-1 + 6^ 
ao = nbl and cq = _ ° (6.7) 

The series for y consists of even terms only. The solution depends upon two 
arbitrary constants only, bo and to, and so is a peculiar solution. 

The system derived for a self-interacting scalar field with an exponential en- 
ergy density evolving in a Robertson- Walker spacetimc containing a separately 
conserved perfect fluid is |g^, ^ 

X = -3x + Ji^ny^ + 7:x [(87 - 2) f^ + 2 (l + 2x2 - ^2^j 

y = -\j\i^xy + iy [(37 - 2) r! + 2 (1 + 2x2 - ^2)] 

(l = -i-in + n [(37 - 2) f7 + 2 (1 + 2a;2 _ y^)] (6.8) 

in a standard notation. This system, with —K in place of k, has already been 
given an extensive treatment in [Q and we simply summarise the results for 
completeness. The leading order exponents are p = q = ^ and r = — 1. Under 
the transformation 

w = 4j;2, v = -2t/ and w = (37 - 2)f} (6.9) 



the dominant terms of (6.8) constitute the system 

u — u{u + V + w) 

i) — V (u + V + w) (6.10) 

w — w {u + V + w) 

which has the leading order behaviour [^ 

u = uqt^ , V ~ voT^ , w — woT~ , (6.11) 

where the coefficients are related according to uo + vo + wq + I = 0. That the 
resonances are r — —1,0(2) confirms that two of the coefficients of the leading 
order terms are arbitrary and these with to provide the three arbitrary constants 
needed for the solution. 

Returning to the original system KM we have that the system passes the 
Painleve test for the weak property. (There is no need to be concerned about 
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incompatibilities at the resonances since the arbitrary constants already appear 
in the leading order terms. There are other patterns of quasisingular behaviour 
similar to those listed in the table in H]. These do not fall within the am- 
bit of the Painleve test and do not affect the possession by ( |6.8| ) of the weak 
Painleve Property and so integrability in terms of functions with only branch 
point singularities (in the case of x and y; il is analytic away from its simple 
pole). 

7 Comments and conclusions 

The analysis of the systems of differential equations in the previous Sections 
enables us to make some more general comments about the Painleve property 
and its possible applications. Firstly, we note that all the systems considered 
here come from reductions of the Einstein equations of general relativity by the 
imposition of certain symmetry and other plausible (from the physical point of 
view) assumptions. That is, they correspond or represent relativistic systems. 
One major characteristic of such systems of differential equations is their covari- 
ance, ie the independence of their properties from the coordinate frame used. 
However, it may be possible that the use of the Painleve analysis for a system 
written in different coordinates leads to different behaviour. Thus we might 
conclude that the Painleve test is not an appropriate one to use when dealing 
with relativistic equations. One may argue then that this test stands at a simi- 
lar level as that of using numerical methods conducive to the determination of 
the so-called Liapunov exponents to decide on the issue of chaoticity of certain 
generic classes of relativistic cosmologies. 

Such an opinion is formulated, as the astute reader will appreciate, somewhat 
tongue in cheek. The literature abounds with studies devoted to the determina- 
tion of closed form solutions of the Einstein equations for various models. The 
variables of the equations are real and, not surprisingly, the solutions sought 
are real. A real solution need not be analytic and yet be quite useful for the 
explication of a specific physical model. Painleve requires solutions which are 
analytic in the complex plane of the independent variable. This is a much 
stronger condition and, as we have remarked above, is superficially in conflict 
with the concept of covariance in the field equations and whatever is derived 
from them. The covariant formulation has no regard for the preservation of 
analyticity. However, the Painleve Property is not preserved under a general 
transformation. Consequently what one should be seeking amongst the plethora 
of possible coordinate representations is that frame in which the Painleve Prop- 
erty can be found, if it is ever to exist for a given system. This, of course, raises 
a new line of research which it is not appropriate to follow in this paper and 
that is to determine which of the possible myriads of coordinate representations 
of the field equations of a given model will be that set of coordinates for which 
the Painleve Property holds. Given that the system is integrable in terms of 
analytic functions in one frame, the lesser demands of Cosmology can be met 
in a wider variety of coordinate systems. 
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